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PROBLEM 


Apply  the  method  of  stationary  phase  to  the  integral  defining 
Woodward's  time-frequency  correlation  function  X(t,  (3)  and  in¬ 
vestigate  the  usefulness  of  a  particular  elaboration  of  this  ap¬ 
proximation  method  as  a  tool  for  attacking  the  signal  design 
problem,  the  latter  being  the  finding  of  such  amplitude  a(t)  and 
phase  4>(t)  modulations  that  a  prescribed  ambiguity  function 
Ix(t,  j3)|2  is  obtained. 


RESULTS 

The  method  is  heuristic  and  is  applicable  to  signals  such  that 
amplitude  a(t)  varies  slowly  relative  to  phase  <J>(t)  over  most  of 
the  signal  duration,  i.  e.  ,  to  phase -modulated  signals.  Those 
regions  of  the  r-P  (range-Doppler)  plane— except  for  a  small 
neighborhood  of  the  origin— on  which  the  ambiguity  function  is 
relatively  large  are  easily  found.  It  is  possible  to  see  that  by 
introducing  finite  discontinuities  in  the  instantaneous  frequency 
d/(t),  some  of  the  ambiguity  volume  (the  total  volume  under 
Ix(t,  p)|2  is  fixed)  can  be  displaced  outward  and  beyond  physi¬ 
cally  interesting  values  of  p.  Further,  the  undisplaced  volume 
can  be  distributed  in  azimuth  about  the  origin.  Hence,  it  seems 
possible  to  approximate  the  ideal  ambiguity  function— sharp  peak 
at  the  origin,  small  values  elsewhere  — in  this  way. 


RECOMMENDATIONS 

The  phase  modulations  arrived  at  in  the  later  sections  of  this 
report,  as  well  as  any  others  obtained  in  the  same  way,  are 
best  regarded  as  constituents  of  candidate  signals.  Additional 
computation— digital,  analog,  or  hand  (using  Erdelyi's  asymptotic 
series)  — is  required  in  order  to  obtain  a  more  accurate  repre¬ 
sentation  of  the  ambiguity  function.  Consideration  of  the  accuracy 
which  may  be  attained,  and  of  the  time  and  effort  required,  indi¬ 
cate  that  the  first  two  are  to  be  preferred  over  the  third. 


1.  INTRODUCTION 


Beginning  with  Woodward  (Ref.  1),  the  time-frequency  correlation 
function 

X(T,/3)  =  J  u(t)u(t  +T)e  **  dt  (1.1) 

has  played  an  important  role  in  the  analysis  of  echo- ranging  systems 
(Ref.  2  to  7).  Much  as  been  written  about  this  function;  Ref.  2  to  4  are 
concerned  with  its  properties,  its  relevance,  and  its  interpretation. 

The  symbols  in  Eq.  1.1  have  the  following  meanings:  the  variable  of 
integration  t  has  the  dimensions  of  time,  t  is  related  to  the  range,  (3 
is  related  to  the  range  rate,  u(t)  is  the  complex  modulation,  i  =  1, 

and  the  overscore  denotes  complex  conjugation.  The  appearance  of  u(t) 
in  Eq.  1.1  results  from  writing  the  physical  transmitted  signal  s(t)  as 

s(t')  =  Re[u(t')e2,lfct  ]  (1.2) 

where  fc  is  the  carrier  frequency  and  t'  is  the  time,  and  from  treating 
the  target  as  a  point  scatterer. 

The  ambiguity  function,  Ix(t,  0)|2,  has  more  physical  significance 
than  X(t,  (3),  and  the  problem  presented  here  is  the  one  of  choosing  the 
modulation  u(t)  so  that  an  acceptable  surface,  Ix(t,  (3) I2  versus  t  and 
p,  is  obtained.  The  qualitatively  ideal  surface  is  usually  taken  to  be  a 
I x (t,  p)|2  which  is  sharply  peaked  at  the  origin  and  relatively  small 
elsewhere.  Aside  from  the  "random"  phase-reversal  signals  (Ref.  3,8) 
and  the  associated  statistical  interpretation  of  Ix(t,  p)l2,  it  appears 
that  a  u(t)  which  yields  an  ideal  Ix(t,  p)  I2  has  yet  to  be  discovered.  As 
might  be  expected,  the  lack  of  such  an  ideal  modulation  has  not  proved 
to  be  crucial.  The  tactical  requirements  of  the  system  may  admit  some 
ambiguity  (Ref.  9)  owing  to  the  presence  of  side  lobes  (additional  peaks 
and  ridges  in  the  Ix(t,  p)|2  surface).  Or  it  may  be  possible  to  resolve 
the  ambiguity  (1)  by  observing  targets  over  a  sufficiently  long  period 
of  time  (Ref.  10)  or  (2)  by  transmitting  more  than  one  kind  of  signal 
(Ref.  4).  On  the  other  hand,  it  may  be  desirable  or  necessary  to  elim¬ 
inate  all  ambiguities  on  each  transmitted  signal.  In  this  last  connec¬ 
tion,  we  note  that,  since  arbitrarily  large  values  of  |t|  and  | (3 1  are  not 
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expected  to  occur,  it  suffices  to  have  lx  (t,  (3)l2  adequately  small,  ex¬ 
cept  for  the  central  peak,  only  on  a  rectangular  domain  of  the  t  -  p 
plane  bounded  by  ±|r|m  and  ±lp|m.  The  particular  elaboration  of  the 
method  of  stationary  phase  that  is  described  below  indicates  how  this 
latter  sort  of  ambiguity  function  might  be  approximated  in  the  case  of 
a  transmitted  signal  consisting  of  a  single  pulse,  as  opposed  to  a  pulse 
burst  or  train  (Ref.  4). 

Although  there  is  a  theorem  due  to  Siebert  and  to  Wilcox  (Ref.  2)  by 
means  of  which  u(t)  can,  in  principle,  be  calculated  if  x(t,  (3)  is  pre¬ 
scribed,  it  appears  that  there  is  neither  a  theorem  leading  directly 
from  Ix(t,  p)|2  to  u(t)  nor  a  theorem  telling  anything  explicit  about  a 
possible  relationship  between  the  magnitude  and  the  phase  of  x(t,  p). 

The  signal  design  problem  must  therefore  be  attacked  by  devious  means. 

Let  us  recall  that  the  method  of  stationary  phase  furnishes  an  ap¬ 
proximate  evaluation  of  the  integral 

F(t')  =J  g(x)e,,,f(X)  dx  (1.3) 

where  x,  f(x),  and  v  are  real,  and  v  is  large  and  positive.  The  rele¬ 
vance  of  the  stationary  phase  approximation  (SPA)  to  the  signal  design 
problem  is  due  to  the  fact  that  the  integral  in  Eq.  1.1,  as  well  as  the 
Fourier  integrals  for  u(t)  and  its  spectrum  U(f),  and  the  integral  ob¬ 
tained  by  applying  Parseval's  formula  to  the  right-hand  member  of 
Eq.  1.1  can  be  cast  into  the  form  of  the  integral  in  Eq.  1.3.  In  the  case 
of  Eq.  1.1,  upon  writing  the  complex  modulation  u(t)  in  terms  of  the 
amplitude  modulation  a(t)  and  the  phase  modulation  <j>(t), 

u(t)  =  a(t)e  (1.4) 

we  obtain 

X(T,  0)  =  f  C  A(t;T)e2l,i*(t’  ’^dt  (1.5) 

where 

A(t;r)  =  a(t)a(t  +  r)  (1.6) 

*(t;T,0)  =  0(t)  -  0(t  +  t)  -  0t  (1.7) 

and  we  have  changed  to  finite  limits  to  agree  with  the  fact  that  we  shall 
be  dealing  with  time-limited  modulations: 
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a(t)  =  0;  t  <  t1(  t  >  tfl 


(1.8) 


The  SPA  also  yields  some  simple,  auxiliary  results  that  are  useful  in 
signal  design.  It  is  in  this  last  respect  that  the  application  of  the 
method  of  stationary  phase  described  and  illustrated  here  differs  from 
previous  applications  (Ref.  2  to  4,  and  11  to  14). 

After  a  brief  review  of  the  SPA  itself  (Section  2),  this  report  de¬ 
velops  (Section  3)  inequalities  involving  t  and  (3  which,  to  within  the 
SPA,  specify  the  domains  of  the  t  -  (3  plane  on  which  Ix(t,  (3)|2  is 
large.1  We  call  such  a  domain  a  side  lobe,  since  the  domain  corre¬ 
sponding  to  the  central  peak  of  Ix(t,  (3)|2  is  excluded  from  these  calcu¬ 
lations.  It  is  then  ascertained  (Section  4),  by  means  of  specific  ex¬ 
amples  and  some  heuristic  considerations,  the  kind  of  phase  modulation 
4>(t)  that  results  in  ( 1 )  the  displacement  of  some  side  lobes  outward  along 
the  (3 -axis  so  that  they  lie  beyond  ±!{3|m  ,  and  (2)  the  distribution  in  azi¬ 
muth  of  those  that  are  contiguous  to  the  origin.  Finally,  the  results 
are  summarized  and  discussed  in  Section  5. 


2.  THE  STATIONARY  PHASE  APPROXIMATION 


The  SPA  is  applicable  to  integrals  such  as  are  given  by  Eq.  1.3 

F(V]  =£  g(x)J't(x)  dx  (2.1/ 

where  f(x)  and  x  are  real  and  v  is  large  and  positive.  This  approxima¬ 
tion  has  been  used  frequently  and  apparently  successfully  in  the  analy¬ 
sis  of  physical  problems.  It  is  discussed  in  many  places,  usually  as 
the  method  of  stationary  phase.  Most  of  the  material  that  follows 
comes  from  the  monographs  by  Copson  (Ref.  15)  and  Erdelyi  (Ref.  16). 
The  basic  idea  is  as  follows:  because  of  the  oscillating  exponential, 
neighboring  contributions  of  the  integrand  tend  to  cancel  except  in  the 
vicinity  of  a  point  Xj  at  which  the  phase  vf(x)  is  stationary  [f'(xj)=  0] 
and  at  the  end  points. 

Copson  assumes  that  f(z)  and  g(z)  are  analytic  functions  of  z  (com¬ 
plex)  in  a  simply  connected,  open  region  D  containing  the  interval  I, 


•The  possibility  of  doing  this  h*s  also  been  noted  by  Vakman  (Ref.  3).  in  connection  with  cubic  phase  modulation 
d-ftloctJ. 
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a  <  x  <  c.  It  then  follows  that  f*(x)  has  a  finite  number  of  zeros  on  I; 
these  are  the  stationary  points  of  the  phase  vf(x).  The  interval  I  is  pre¬ 
sumed  to  be  divided  into  a  finite  number  of  subintervals,  in  each  of 
which  f’(x)  does  not  vanish,  or  f'(x)  vanishes  only  at  the  left-hand  end 
point,  or  f'(x)  vanishes  only  at  the  right-hand  end  point.  These  three 
cases  are  treated  separately,  and  the  results  are  as  follows: 


If  f(x)  has  no  stationary  point  in  a  <  x  <  K,  then 

»*<x)  ,  _  g(*Q .  ivf<*)  g(tt)  i»f(a)  /J_\ 

“  '  iUf'(OL)e  +  °U  2/ 


.  rK  ,  ,  iv^x)  ,  g(*0 

G  (K)  g(x)e  dx  =  e 


(2.?) 


as  v  -*■  oo. 


If  f(x)  has  one  stationary  point  ina<x</c,  namely  at  x  =  a,  then 

_  ,  ,  f*  .  .  ivf(K) 

G  (v)  =  I  g(x)e 

J  U 

g(a)e  +°W„  (  ,  ) 

as  v  ■*  »,  where  the  upper  (lower)  signs  apply  when  f"(a)  is  positive 
(negative). 

If  f(x)  has  one  stationary  point  in  a  <  x  <  k  ,  namely  at  x  =  k,  then 
G  (V)  =Ja  g(x)e  dx 

=  «2-4» 

as  v  -*•  oo  ,  where  the  upper  (lower)  signs  apply  when  f"(K)  is  positive 
(negative). 

Although  these  results  will  suffice  for  most  of  the  discussion  in  the 
next  two  sections,  a  more  general  theorem  is  useful  (1)  for  those  cases 
where  either  f"(a)  =  0  or  f’(»c)  =  0;  (2)  for  establishing  that  the  right- 
hand  members  of  Eq.  2.2  through  2.4  are,  in  fact,  the  dominant  terms 
of  an  asymptotic  expansion;  and  (3)  for  obtaining  more  accurate  asymp¬ 
totic  approximations.  This  theorem  is  derived  by  Erdelyi  and  includes 
stationary  points  of  a  more  general  character.  It  is  assumed  that  the 
number  of  these  more  general  stationary  points  in  the  interval  I  is 
finite,  and  that  I  is  broken  up  as  before  except  that,  for  definiteness, 
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f(x)  is  monotonic  increasing2  on  a  <  x  <  k.  The  theorem  reads  as 
follows: 


If  \  >  0,  |i  >  1;  g(x)  is  N  times  continuously  differentiable  for 
a  <  x  <  k;  f(x)  is  differentiable  and 

f'(x)  =  (x  -  a)fi  1  (  k  -  x)a_1  fu'x)  (2.5) 


where  p,  a  >  1,  and  fi(x)  is  positive  and  N  times  continuously  differ¬ 
entiable  for  a  <  x  <  k;  then 


f  g(x)ell,f<X)  (x  -  a)X  *  (k  -  x)^  1  dx  =  B(l/)  -  A(v) 
J  a 

where 


and 


A(V)  ~ 


iri(n+X)/2p  ~(n+X)/p  ivf(a) 
e  v  e 


B(V) 


N-l 

-  z 


n=0 


l(n)(0)  r 

n!  a 


-iri(n -{n +p)/o 
e  v  e 


as  v  -*•  oo,  with  k(u)  given  by 

k(u)  =  gi(x.)u'k  ^ 
up  =  f (x)  -  {(a) 

gi(x)  =  g(x)(x  -  a)*"1  (K  -  x)u  1 


(2.6) 


(2.7) 


(2.8) 


(2.9) 


and  t  (v)  given  by 


The  notation 


*(v)  =  giWv1-*1  ~ 
v°  =  f(K )  -  f(x) 


N-l 

A{u)  ~  £  a„0n(l/)  as  V  -  oo 
n«0 


(2.10) 


^  If  i(x)  is  monotonic  decreasing,  replace  i  with  -i  and  f(x)  with  -f(x). 
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means  that 


N-l 

A(f)  =  £  antf’nU')  +  °(0N-1 )  as  v  -  00 

r.=0 

The  presence  or  absence  of  a  stationary  point  within  the  interval  of 
integration  is  manifested  in  the  values  of  p  and  cr.  If  there  is  no  sta¬ 
tionary  point,  p  =  <r  =  1;  otherwise,  either  p  or  cr,  or  both,  is  greater 
than  unity. 


3.  GENERAL  CONSIDERATIONS 


The  use  of  the  SPA  on  Eq.  1.5  is  now  considered  in  more  detail.  No 
loss  in  generality  results  from  choosing  the  time  origin  so  that 


T  T 

ti  -  -  2  ,  ts  =  +  2 

(3.1) 

We  then  have 

a(t)  =  0  for  t<  -p  t>j 

(3.2) 

from  Eq.  1.8. 

Also,  because  of  the  symmetry  condition, 

1  X(-T.  -  0)|  =  |x(T,0)| 

(3.3) 

we  may  take 

T  >  0 

(3.4) 

throughout.  Equation  1.5  then  becomes 

X<T,0)  A(t;T)e2"‘,,:''«dt 

J-T/2 

(3.5) 

where 

*(t;T,p)  =  0 (t )  -  0 (t  +  T)  -  (3t 

(3.6) 

and 

A(t;T )  =  a(t)a(t  +  t) 

(3.7) 

We  note  that 

X(t,  0)  =  0  for  r  £  T 

(3.8) 

The  simplest  case  obtains  when  <Jj(z)  and  a(z)  are  analytic  functions 
of  z(Rez  =  t)  on  a  simply  connected,  open  region  D  containing  the  t- 
interval  [-T/2,  T/2].  Since  this  assumption  conflicts  with  Eq.  3.2,  we 
have  to  suppose,  for  the  purpose  of  applying  the  first  three  theorems 
given  in  the  previous  section,  that  a(t)  is  analytically  continued  to  the 
left  beyond  -T/2  and  to  the  right  beyond  +T/2.  The  theorems  of  the 
previous  section  then  tell  us  that  the  SPA  will  be  a  good  approximation 
if,  aside  from  stationary  points,  the  phase  2tT4>(t;  t,  (3)  varies  rapidly 
with  t  compared  to  A(t;  t)  on  [-T/2,  T/2].  We  assume  this  to  be  true 
henceforth,  and  it  then  follows  from  Eq.  3.6  and  3.7  that  the  phase 
modulation  4>(t)  must  vary  relative  to  the  amplitude  modulation  a(t)  on 
[-T/2,  T/2].  In  a  specific  case,  a  factor  v  can  be  taken  out  of 
2ir4>(t;  t,  (3).  These  theorems  also  tell  us  that  the  ambiguity  function 
|x(t,  P)l2  will  be  larger  by  a  factor  v  when  the  stationary  point(s)  lie 
within  the  interval  of  integration;  i.  e.  ,  when 


T 

2 


£ 


with  the  stationary  point  t.j  satisfying 


(3.9) 


p*(t:p  g)]  =  jV(t)  .  0'(t  +  T)  -  fi]  =o  (3.10) 

L  ot  J  t=tj  L  Jt=tj 

Hence,  tj  is  a  function  of  t  and  p,  and  Eq.  3.9  may  be  written  as 

-j*t  ,M)s|-r  (3.11) 

The  equations 


(r,0)  =  -j 

t)(T,P)  =J  -  T 


(3.12) 


define  the  loci,  p  versus  t,  cori  esponding  to  a  stationary  point  at  the 
ends  of  the  interval  of  integration  [-T/2,  (T/2)  -  t].  According  to 
Eq.  3.10,  these  loci  are  given  by 


P  =  M*)  "  )]j=_T/2  H  fl(T) 


(3.13) 


and 


P  =  [^'(t)  -  *'(t  +  T)]t=T/2_T  =  f s(T) 


(3.14) 
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We  shall  refer  to  fj  (t)  and  f2(*r)  as  the  terminal  curves.  When  the 
stationary  point  tj  is  within  the  interval  of  integration,  the  loci  are 
given  by 

0  =  [0'(t)  -  0'(t  +  T)]  ,  -7<tj  <f-T  (3.15) 

from  Eq.  3.10.  It  is  obvious  that  all  these  loci3  pass  through  the  origin, 
t  =  p  =  0.  They  also  intersect  at  t  =  T,  as  may  be  seen  by  observing 
that  there  is  then  only  one  possible  value  for  tj,  namely  tj  =  -T/2. 
Collectively,  these  loci  constitute  the  side  lobe,  the  boundaries  of 
which  will  be  either  certain  members  of  the  family  of  loci  which  en¬ 
close  all  the  others  or  they  will  be  envelope  curves. 


In  this  connection,  we  consider  next  the  differential 


d0  =  [0"(t)  -  0"(t  +  r)],=tj  dtj  = 


(3.16) 


which,  for  fixed  t,  gives  the  displacement  dp  involved  in  moving  from 
one  locus  to  another  which  is  close  by.  By  looking  at  the  sign  of 
[324>(t;T,  p )/ 9t 2  ]t=tj  as  the  stationary  point  tj  moves  into  the  interval  of 
integration,  we  shall  be  able  to  see  how  the  loci  (Eq.  3.15)  are  dis¬ 
posed  relative  to  the  terminal  curves  f  j  (t)  and  f2(T).  Clearly,  the 
terminal  curves  and  the  comments  relative  to  Eq.  3.15  apply  to  any 
other  stationary  points  which  may  exist  in  addition  to  tj.  Hence,  the 
extent  of  the  side  lobe  may  be  determined  in  this  way  regardless  of 
how  many  stationary  points  are  present.  However,  the  existence  of 
more  than  one  stationary  point  may  influence  the  magnitude  of  x(t,  p) 
on  the  side  lobe:  by  applying  Eq.  2.3  and  2.4  to  Eq.  3.5,  we  obtain 


Xlr.0)  -  £  '3.17) 


with  tj  =  tj ( t,  P)  a  solution  of  Eq.  3.10.  It  also  has  the  effect  of  making 
the  SPA  less  accurate  as  t  -*  T  from  below  and  the  multiple  stationary 
points  coalesce  (a  circumstance  which  is  not  considered  important,  in 
view  of  Eq.  3.8). 


The  simplest  situation  obtains  when  9 2 <t>( t;  t,  p ) / 9 1 2  is  independent 
of  t  and  has  the  same  sign  on  0  <  t  <  T;  say 


;t,/3) 

dt2 


0  <  T  <  T 


(3.18) 


3Thc  fact  that  a  parametric  family  of  curves  is  generated  by  varying  the  stationary  point  was  noted  by  Sorkin 
(Ref.  13)  but  not  developed  in  detail. 
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(We  note  in  passing  that  9 2 4>( t;  t,  (3)/9t2  is  always  independent  of  (3.  )  If, 
in  addition, 


fi(T)  <  f2(r),  0  <  r  <  T  (3.19) 

then  the  side  lobe  is  specified  by 

fx (r )  s  0  £  f2(r),  0  <  r  <  T  (3.20) 

an  inequality  which  is  equivalent  to  Eq.  3.9  and  which  occut:  in  Ex¬ 
ample  1  in  the  next  section. 

The  question  arises  whether  the  inequality  sign  could  be  reversed 
in  either  (3.18)  or  ( 3. 1 9 )»  for  in  that  event  there  would  be  no  side  lobe 
at  all,  a  desirable  state  of  affairs.  However,  this  would  lead  to  the 
contradiction  that  fi(T)  and  f2(T)  are  not  members  of  the  same  family  of 
curves,  there  being  no  way  of  getting  from  one  to  the  other  by  varying 
tj.  The  possibility  that  one  of  the  signs  may  be  reversed  over  just 
0  <  t  <  T/2,  for  example,  may  be  disposed  of  on  the  same  grounds. 

The  more  complicated  cases  arise  when  92 «t»( t;  t,  (3 ) / 9t2  depends 
upon  both  t  and  r.  An  enumeration  and  analysis  of  these  has  not  been 
attempted  here  since  the  general  procedure  seems  clear.  Starting 
with  f i ( t)  (the  member  of  the  family  of  loci  corresponding  to  tj  =  -T/2), 
the  evolution  of  the  remaining  members  is  ascertained  by  repeated  ap¬ 
plication  of  Eq.  3.16  as  tj  moves  from  -T/2  to  T/2  -  t.  If  more  con¬ 
venient,  it  is  possible  to  start  with  f2(T),  in  which  case  tj  evidently 
moves  from  (T/2)  -  t  to  -T/2,  and  dtj  in  Eq.  3.16  is  negative.  It  is 
not  uncommon  for  92<l>(t;  t,  P)/9t2  to  change  sign,  perhaps  at 
-T/2<$<T/2-t.  If  Eq.  3.19  holds,  and  the  change  in  sign  is  from 
plus  to  minus,  then  the  side  lobe  will  be  given  by 

fx(T)  s  |8  s  f(T),  0  <  t  <  T  (3.21) 


where 


f(T)  =  [0'(t)  -  0'(t  +T)]t.} 

is  the  equation  for  the  envelope  curve  (Ref.  17),  Examples  2  and  3 
below  are  instances  of  this  sort. 

Suppose  now,  contrary  to  the  assumption  following  Eq.  3  8,  that 
the  amplitude  a(t)  has  finite  discontinuities  at  t  =  ±T/2.  Such  a  case 
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is  of  interest  as  a  mathematical  idealization  of  certain  physical  wave¬ 
forms.  Let  h(t)  be  defined  by 

h(t)  =  a(t)  for  1 1 1  <  -- 

T 

h(t)  =  analytic  continuation  of  a(t)  for  |t|  ^  —  (3.22) 

Then,  by  using 

A(t;  T)  =  h(t)h(t  +  T)  (3.23) 

in  place  of  Eq.  3.7,  the  previous  results  are  again  obtained. 

When  there  are  finite  discontinuities  in  a'(t),  4>(t),  <J>'(t),  etc.  ,  the 
interval  of  integration  [  — T/2,  (T/2)  -  t]  is  partitioned  so  that  the 
points  of  discontinuity  coincide  with  the  limits  of  the  integrals  making 
up  Eq.  3.5,  and  then  analytic  cc  tinuations  of  a(t)  and  4>(t)  are  intro¬ 
duced  in  order  that  the  theorems  of  Section  2  can  be  applied  to  each  of 
the  constituent  integrals.  Of  course,  the  analytically  continued  func¬ 
tions  have  a  strictly  formal  role;  all  that  needs  to  be  done  in  a  calcula¬ 
tion  is  to  avoid  integrating  through  a  discontinuity.  Results  similar 
to  Eq.  3.11,  3.13,  3.14,  3.15,  and  3. 1 6  then  apply  for  each  segment  of 
the  original  interval  [-T/2,  (T/2)  -  t]. 

In  the  case  of  Erdelyi's  theorem  (Eq.  2.6  through  2.8),  neither 
a(t)  nor  <j>(t,  need  be  analytic.  It  has  not  been  possible  to  see  whether 
this  weaker  hypothesis  has  any  desirable  consequences  insofar  as  the 
signal  design  problem  is  concerned,  except  that  the  theorem  gives  a 
finite  result  in  some  cases  in  which  Eq.  2.3  and  2.4  do  not.  In  particu¬ 
lar,  on  the  (3-axis,  82<|>(t;  t,  0)/8t2  vanishes  identically,  and 

X(T,/3) 

0-0 

according  to  Eq.  2.3  and  2.4,  while  Erdelyi's  theorem  gives  a  finite 
result,  which,  for  the  case  of  a  rectangular  a(t),  turns  out  to  be  exact. 
This  circumstance  is  useful  in  a  supplementary  way  in  that  degenerate 
side  lobes  can  be  defined  (Eq.  4.6  and  4.7  below). 

The  SPA  fails  completely  at  the  origin  because  4>(t;  0,  0)  =  0.  But 
it  is  known  from  other  considerations  that  lx  I2  attains  a  maximum 
there.  Also,  its  behavior  in  a  small  neighborhood  of  the  origin  can 
be  studied  by  means  of  a  Taylor  series  expansion. 


Thus,  our  use  of  the  SPA,  principally  to  map  out  those  regions  of 
the  plane  on  which  lx(T*  P)l2  is  large,  excludes  the  central  maximum 
and  a  narrow  strip  about  the  |3-axis.  Along  the  (3-axis,  one  can  use 
either  Erdelyi's  theorem  or  the  exact  equations 

X(T .  /?)  =  /_  [a(t)f  e’2'^dt=J“  A(f)  A(f+/J)df  (3.24) 

which  follow  from  Eq.  3.5  and  Parseval's  formula.  Here,  A(f)  is  the 
spectrum  of  a(t). 

It  is  natural  to  inquire  into  the  possibility  of  attacking  the  signal 
design  problem  synthetically  within  the  present  context,  i.  e. ,  by 
specifying  the  bounding  curves  of  the  side  lobes,  then  working  toward 
the  modulation,  Eq.  1.4.  It  is  easy  to  convince  oneself  that  this  pro¬ 
cedure  is  ambiguous.  In  particular,  the  integrations  with  respect  to  t 
that  are  necessary  lead  to  functions  of  t  which  may  be  specified  arbi¬ 
trarily.  Hence,  it  seems  more  satisfactory  to  work  toward  the  place¬ 
ment  of  the  side  lobes  by  trial  and  error. 


4.  SOME  EXAMPLES 


In  this  section,  the  elaboration  of  the  SPA  described  above  is  ap¬ 
plied  to  several  different  phase  modulations  (j>(t).  The  amplitude 
modulation  a(t)  is  unspecified,  except  that  it  is  assumed  to  be  such 
that  the  theorems  of  Section  2  are  applicable,  and  to  vary  slowly  com¬ 
pared  to  6(t)  over  the  duration  T  of  the  signal.  The  supposition  that 
Ix(t,  (3)  |2  is  large  on  the  side  lobes  (stationary  points  within  the  in¬ 
terval  of  integration)  and  small  elsewhere  (except  for  the  central  lobe) 
will  then  be  admissible.  For  convenience,  we  introduce  the  function 
6(t;  t): 

®(t;r)  =  0'(t)  -  0'(t  +  T)  (4.1) 

The  two  common  modulations, 

0 (t )  =  c  (monotone)  (4.2) 

and 

0 ( t )  =  j  csts  (linear  FM4)  (4.3) 


*  Frequency  modulation. 
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where  c  and  c2  are  constants,  lead,  respectively,  to  the  equations 


and 


It  , 

ir^=  -caT  -p  = 0 


(4.4) 

(4.5) 


for  the  stationary  point,  Since  t  does  not  appear  in  these  equations, 
there  are  no  stationary  points  of  the  ordinary  kind,  and  Eq.  2.3  and  2.4 
do  not  apply.  On  the  other  hand,  Erdelyi's  theorem  yields  a  finite  re¬ 
sult,  which,  it  turns  out,  is  exact  for  a  rectangular  a(t).  Now,  in  the 
case  of  a  rectangular  a(t),  Ix(t,  P)l2  decreases  in  directions  orthogonal 
to  the  lines 


and 


(3  =  0  (monotone) 

(4.6) 

fi  -  -csT  (linear  FM) 

(4.7) 

It  is,  therefore,  not  unreasonable  to  consider  that  these  equations  con¬ 
stitute  a  degenerate  sort  of  side  lobe. 


Example  1. 

A  slightly  more  complicated  phase  modulation  is 

0(t)  -  \  c2t3  +  }  c3t3  (4.8) 

which  may  be  called  linear  plus  quadratic  FM.  The  term  linear  in  t  has 
been  omitted  because  it  has  no  effect  on  Ix(t,  (3)|2.  We  have 

*(t ;r,fi)  =  [j  cst3  +  y  c3  t3] 

-  [y  cs(t  +  T)3  +  j  c3(t  +  r)3]  -  fit  (4.9) 

=  -[c3t3  +  (cs  +  Zc3t)  T  +  fi]  (4.10) 


®(t;r)  =  -[c3r3  +  (c8  +  2c3t)  t] 


(4.11) 


and 


a24>(t;r,  8) 

dt3 


=  -2c3r 


(4.12) 


i 

! 


i 


t 


i 


2 


F  rom  the  equation 


d$(t;T,  8) 


at 


=  0 


we  find  that  there  is  just  one  stationary  point  tj,  given  by 

tj  =  (c3T2  +  c2T  +  0) 

I  he  terminal  curves  f ,  (t)  and  f2(T)  are  given 


by 


and 


fi(r)  =  [0(t,;r)]  T/2 

=  -c*[T+2  (^-T)]S+^(fl-T)' 

fe(r)  =  [®(t,;r)lt^T/1_T 

*«•  [T-1(^+T)]S-  4 


(4.13) 


(4.14) 


=  C4T-2  (t  +  T)J  -“4  (t+T)  (415) 

and  the  remaining  members  of  the  family  constituting  the  side  lobe  by 
0  =  -c3r3  -  (Cjg  +  2c3tj)r,  -  \<  tj  <  \  -  r  (4.16) 

Suppose,  for  definiteness,  that 


c3<  0,  &  -  T  >  0  (4.17) 

c3 

The  loci  of  the  terminal  curves  then  are  as  sketched  in  Fig.  1.  On  the 
interval  0  <  t  <  T,  fj(T)  <  f2(T),  and  from  Eq.  4.12  and  4.17, 

824>(t;  t,  (3)/3t2  >  0.  If  Eq.  3.16  is  now  presumed  to  be  applied  repe- 
titiously,  starting  from  tj  =  -T/2  and  the  curve  f|(T),  we  can  see  that 
the  inequality  (3.20)  applies,  and  that  the  side  lobe  of  !x(t,  0)l2  for 
t  >  0  is  as  indicated  by  the  shaded  region.  The  loci  (4.16)  lie  within 
this  region  on  the  interval  0  <  t  <  T. 


For  convenience,  such  a  diagram  is  called  an  ambiguity  diagram 
here.  The  part  for  t  <  0  is  obtained  by  reflecting  the  shaded  region  in 
the  origin  (see  Eq.  3.3).  The  central  lobe  is  not  shown. 

Example  2. 

In  this  example,  we  take 

0>(t)  =  ^  ct4  (4.18) 
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FIG.  1.  Ambiguity  Diagram  for  Linear  PIu>  Quadratic  FM,  Eq.  4.8;  r>0. 
|x(r,  fi)P  it  large  on  the  thaded  region. 


which  may  be  called  cubic  FM,  and  for  which  the  following  results  are 
obtained: 

*<t;T,0)  =^ct4  -  |c(t  +t)4  -  fit  (4.19) 


=  -cr(3t2  +  3tr  +  r2)  -  0 


(4.20) 


@(t  ;t)  =  -crOt2  +  3tT  +  t2) 


(4.21) 


and 


324>(t;r(l 

dt2 


-3cr(2t  +  r) 


From  Eq.  4.20  it  is  seen  that  there  are  two  stationary  points, 


(4.22) 


(4.23) 


The  terminal  curves  are 

f,(r)  =  [9<t,;T)]v  T/2  =  -ct[(t  -  ”)  +3(1)]  (4-24> 

£s(t)  =  =  -cr[(r  -  ”)  +  3  (7)  ]  (4  25) 


14 


Since  the  loci  of  f,  (t)  and  f 2(t)  coincide  and  82 4»( t;  t,  (3)/9t2  depends 
upon  t,  the  ambiguity  diagram  for  this  example  will  have  a  character 
entirely  different  from  that  of  the  previous  case.  For  use  in  Eq.  3.16, 


f  S24>(t ,;t,  g)1 
L  3tj2  J  tj«-T/2 


3cr(T  -  r) 


(4.26) 


Taking 


[ 


afMiiL-Jil 

dtj2  Jtj-T/2-r 


-3cr(T  -  T) 


c  >  0 


(4.27) 

(4.28) 


and  referring  to  Eq.  3.16  and  4.26,  it  is  found  that,  as  a  stationary 
point  moves  from  -T/2  into  the  interval  [-T/2,  (T/2)  -  t],  the  loci 

/3  =  [®(t;r)]t=t_  =  -cT<3t,2  +  3t,7  +72),  -\<  t ,<y-r  (4.29) 

are  displaced  upward  from  fj(T).  Similarly,  Eq.  3.16  and  4.27  reveal 
that,  as  a  stationary  point  moves  from  (T/2)  -  r  into  the  interval 
[-T/2,  (T/2)  -  t],  the  loci  (4.29)  are,  again,  displaced  upward  from 
f  l  ( t )[  =A f 2 (T)]<  It  is  therefore  necessary  to  determine  the  envelope 
curve  f  (t )  corresponding  to  the  solution  t  of 


524»(t;r,  B )  _ 

at2 


-3cr(2t  +  r)  =  0 


(4.30) 


By  inspection,  this  solution  is 


U  -I 


(4.31) 


from  which 


f  (r)  =  [®  (t;T )]  a  =-tt3 


't=t 


(4.32) 


and 


p^(t;r,fl)j 


A 
t  =  t 


=  - 6cr  <  0 


It  follows  that  the  side  lobe  is  specified  by 


(4.33) 


fi(7)  s  0  s  f(T),  0  <  7  <  T 


(4.34) 


The  ambiguity  diagram  for  t  >  0  and  with  the  central  lobe  excluded 
appears  in  Fig.  2. 


FIG.  2.  Ambiguity  Diagram  for  the  Modulation  of  Eq.  4.18,  Cubic  FM. 


It  may  be  worth  noting  that  explicit  use  was  not  made  of  the  expres¬ 
sion  for  the  stationary  points  given  in  Eq.  4.23.  That  equation  is  in¬ 
cluded  only  for  completeness;  were  Ix(t,  (3)  I2  to  be  calculated  on  the 
side  lobe,  it  would  be  needed  for  substitution  into  Eq.  3.17.  Similar 
comments  apply  to  Eq.  4.13  in  the  preceding  example. 

Example  3. 

Because  the  case  of  sinusoidal  FM 

0(t)  =  cx  sin  Cgt 

is  somewhat  awkward  to  discuss  for  arbitrary  values  of  c2 ,  a  specific 
and  convenient  value  is  assumed: 

0(t)  =  cx  sin  2ff  y  (4.3E 

From  this  there  follows 


t  t  +  T 

4>(t;r,/3)  =  cx  sin  2n  —  -  cx  sin  2v  -  8 1 


t+; 


d4>(t;r,i3)  4jfc,  .  t  .  _  2 

- - -  =  ^  1  sin  if  —  sin  2lf  — — —  -  a 

at  T  T  T  H 


(4.36 

(4. 31 


T 

© (t ;T )  =  4*C*  sin  j r  zz  sin  2rr  — ~~ 


(4.38 


fi(T)  =  [®(tJ;r)Jt  =  T/2 


2  ffC! 

T 


[l  +  COS  Zv{l-  j)] 


(4.39) 


f=(r)  =  [0(tJ;r)]t^T/2_T  =  -fx(r) 

d24>(t,r,3)  ^/tt\s  .  t 

- gp -  =  8  cx  sin  TT  —  cos  2^ 


‘  +  l 


(4.40) 

(4.41) 


and 


rd2*(t,;T,g)] 

L  dt2  Jy-T/2 


p2<j>(t,;T,fl)l 

L  atj2  JtpT/2-r 


/2tr \ 3  .  -  T 

-Ci\—  /  sm2ir  — 

p2$(t,;T,gn 

L  Stf  Jtp-T/2 


(4.42) 


(4.43) 


Since  the  equation  for  the  stationary  points,  94>(t;  t,  P)/9t  =  0,  is 
transcendental,  the  evaluation  of  the  ambiguity  function  would  be  diffi¬ 
cult  in  this  case.  The  ambiguity  diagram  can,  of  course,  be  deter¬ 
mined  with  less  effort. 


From  Eq.  4.39,  4.40,  4.42,  4.43,  and  3.16,  we  can  construct  the 
preliminary  diagram  shown  in  Fig.  3  where,  for  definiteness,  we  have 
taken 


cx  <  0  (4.44) 

The  arrows  in  Fig.  3  indicate  the  directions  in  which  neighboring 
members  of  the  family  of  loci  making  up  the  side  lobe  are  to  be  found. 
The  equation  for  these  loci  is 

T  'T’  'T’ 

P  =  [®(t;r)]t=t_  =  — p1  sin  v  —  sin  2tt  — — ,  -  —  <  ts  <  —  -  T  (4.45) 

The  points  Pj  and  P2  seem  to  have  a  sort  of  singular  character,  and 
we  put  t  =  T/2  in  Eq.  4.45  to  find  ort  what  happens  there.  This  gives 

P  =^cos  2tt  -p<tj<0  (4.46) 


which  tells  us  that  the  point  Pj  moves  downward  to  P2  as  the  parameter 
tj  moves  through  its  allowed  range. 
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FIG.  3.  Displacements  From  the  Terminal  Curves  Corresponding  to  Small 
Increments  in  the  Parameter  tj  (Stationary  Point);  Sinusoidal  FM  of  Eq.  4.35. 


In  order  to  elucidate  matters  further,  we  look  for  the  values  t  at 
which  824»(t;  t,  P)/9t2  vanishes.  By  so  doing,  the  envelopes  of  the  loci 
(4.45)  can  be  determined.  It  is  sufficient  to  look  at  the  time-dependent 
factor  in  Eq.  4.41.  This  vanishes  when 


t  =  ~  £(2m  +  1)  -  rj  ,  m  =  integer  or  zero 

By  combining  this  with 

then  using  0  <  t  <  T  (Eq.  3.4  and  3.8),  we  arrive  at 

1 2m  +  1  |  <  2  (l  -  m  =  0,  ±1,  ±2,  ••• 


(4.47) 


(4.48) 


This  inequality  has  the  solutions  m  =  0  or  -1  for  0  <  t  <  T/2,  and  no 
solutions  for  T/2  <  t  <  T.  Upon  using  these  results  in  Eq.  4.47,  then 
substituting  into  Eq.  4.41,  it  is  found  that 

d2»(t;r,l)  =  0att  =  ±?-T.  0  <  T  <  T 


at2 


(4.49) 


I 


) 

! 


I 


I 
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Hence,  for  T/2<t<  T  the  sense  of  the  displacements  shown  in  Fig.  3  per¬ 
sists  throughout  the  allowed  interval  for  tj.  For  0<t<  T/2,  the  envelope 
curves  are  found  by  substituting  t  =  ±T/4  —  t/2  for  tj  in  Eq.  4,45: 


T 

sm  V  — 


&(T)  =  -fi(T) 


(4.50) 


The  sense  of  the  displacements  for  0  <  t  <  T/2  shown  in  Fig.  3  reverses 
at  these  curves. 


The  results  may  now  be  combined  so  as  to  obtain  the  inequalities 
specifying  the  side  lobe, 

fa(T)  £  0  *  fi(T),  0  <  T  £  ~ 

fg(r)  £  P  £  Mr),  -|£T  £  T  (4.51) 

and  the  corresponding  ambiguity  diagram,  shown  in  Fig.  4.  As  before, 
this  diagram  is  for  r  >  0;  the  complete  diagram,  except  for  the  central 
lobe,  may  be  obtained  by  inversion  in  the  origin. 


FIG.  4.  Ambiguity  Diagram  for  the  Sinusoidal  FM  of  Eq.  4.35. 


In  the  examp  considered  thus  far,  the  phase  modulation  <|>(t) 
satisfied  the  hype  eses  of  the  theorems  stated  in  Section  2,  i.  e.  ,  4>(z), 
with  Rez  =  t,  was  analytic  on  an  open,  simply  connected  region 
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containing  the  t-interval  [  -T/2,  T/ 2] .  In  each  case,  it  was  found  that 
the  ambiguity  function  Ix{t,  |3)|2  had  side  lobes  connected  to  the  origin 
and  extending  outward  to  t  =  ±T.  The  discussion  in  Section  3  following 
the  inequality  (3.20)  indicates  that  this  sort  of  result  will  always  be 
obtained. 

We  are  thus  led  to  consider  modulations  <j>(t)  which  are  not  analytic 
as  a  possible  means  to  obtain  an  ambiguity  function  all  or  part  of  whose 
side  lobes  are  displaced  outward  beyond  the  domain  of  values  of  t  and 
P  which  are  expected  to  occur  in  practice.  As  pointed  out  in  Section  1, 
such  a  Ix(t,  P)l2  would  be  a  desirable  approximation  to  the  ideal. 

Example  4. 

To  this  end,  let  us  consider  the  modulations  given  by 

a(t)  =0,  t  <  0,  t  >  T  (4.52) 


<t>,(t). 

°**<F 

d,(t)  = 

(4.53) 

4>j(t), 

Z<  t£  T 

2 

<j)j(t)  =  C  |Q 

*■  cnt  + 

1  ,  1  3 
-  c12t  +  3  C13t 

•M*)  =  C20 

f  C2jt  + 

1  2  1  3 

2  c22t  +  3  c23t 

(4.54) 

The  origin  of  time  has  been  shifted  because  a  certain  unifoimitv  in 
some  of  the  expressions  which  follow  is  thereby  attained.  We  con¬ 
tinue  to  take  t  >  0.  The  sketch  in  Fig.  5  may  be  helpful  in  determining 
the  subdivisions  of  the  interval  of  integration  [0,  T  -  t]  .  Referring 
also  to  Eq.  1.5,  we  have  for  0  <  t  <  T/2, 


x<r.m=/0 


ran-. 


Ait  ;r  ie2l,i*,l<1;T’^  dt 


[Tl2  . . '’dt.  fT" All;rleI'l*3!"’’-'ldt 

J  T/2 


+ 


(4.55) 


FIG.  5.  Aid  for  Determination  of  Limits  in 
Eq.  4.55. 


and  for  T/2  <  t  <  T, 


f T  ^  ,  2* i#12(t;r,lO  . 

X(t,0)=|  A(t  ;r)e  dt 

0 


Here,  the  functions  <|>„ ,  etc.,  are  given  by 
<*>u(t;T ,  fi)  =  0i (t)  -  0i (t  +T)  -  fit 

=  -[cur  +  \  C„(2tr  +  Ts>  +  }  c13(3t»T  +  StT3  t  Ta>] 

4»ia(t;r.  /3)  =  0x(t)  -  0a(t  +  t)  -  fit  =  (c10  -  c30)  +  (cu  “ 

+  J  (CIS  -  C8s)tS  +  J  (c13  -  C23)t3 
-  [cslr  +  \  Cgs(2tr  +r3)  +|  ca3(3rt  +  3tr=  +  T3)_ 


<S»aa(t;r ,  ft)  =  03(t)  -  0a (t  +  T)  -  fit 


=  -[caiT  +  |  Caa(2tT  +  T2)  +  |  c*3(3t?T  +  3tTa  +  T3)]  -  fit  (4.59) 


and,  as  before. 


A(t;T)  =  a(t)a(t  +  T) 


(4.60) 


Were  the  calculation  of  Ix(t,  P)l2  to  be  done  completely,  it  would  be 
necessary  to  take  their  relative  phases  into  account  when  adding  up  the 
three  terms  of  Eq.  4.55.  This  raises  the  possibility  of  contriving  a 
destructive  interference.  However,  such  an  endeavor  would  succeed 
only  if  the  constituent  side  lobes  were  coincident,  which  is  unlikely, 
and  this  course  will  not  be  pursued  here.  Since  the  functions  «J>n(t;  t,  (3) 
and  4>n(t;  t,  p)  have  the  same  character  as  the  function 4>(t;  r,  p)  pre¬ 
viously  dealt  with,  the  side  lobes  contributed  by  the  first  and  third 
integrals  in  Eq.  4.55  will  be  similar  to  those  encountered  before,  ex¬ 
cept  that  they  will  extend  outward  from  the  origin  only  to  t  =  T/2.  The 
second  integral  in  Eq.  4.55,  as  well  as  the  integral  in  Eq.  4.56,  offers 
the  possibility  of  yielding  a  different  kind  of  result.  For  each  of  these 
integrals  we  have,  from  Eq.  4.58, 

34>ia(t;T,  fi)  r,  -  -  -21 

-  -  UC11  -  C21  -  C22T  ■  °S3T  J 

+  (c13  -  Caa  -  2 c33t )t  +  (c13  -  ca3)t2]  -  fi  (4.61) 


=  (ci2  '  caa  ~  c23^ )  +  2(c13  -  ca3)t 


(4.62) 


Starting  with  these  two  equations  and  proceeding  as  in  the  previous  ex¬ 
amples,  we  want  to  see  whether  the  constants  Cjj  can  be  chosen  in  such 
a  way  that  the  side  lobe  corresponding  to  Eq.  4.56  is  displaced  outward 
along  the  p-axis  beyond  the  maximum  expected  values  ±lplm  .  If  it  is 
possible  to  do  this,  then,  to  within  the  SPA,  the  t  -  p  plane  will  be 
effectively  free  of  ambiguity  for  t  >  T/2.  In  view  of  the  remarks  above 
concerning  interference  among  the  terms  of  Eq.  4.55,  we  should  simi¬ 
larly  try  to  displace  the  side  lobe  corresponding  to  the  second  term 
there. 


fV2  A(t;T)e2wi*u<t:r,<,)dt 
J(T/2)-t 


(4.63) 
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This  displacement  of  side  lobes  involves  an  additional  benefit  which  is 
due  to  the  fact  that,  to  within  the  narrow- band  approximation,  the 
volume  under  the  Ix(t,  P)l2  surface  is  independent  of  the  waveform  and 
depends  only  upon  the  signal  "energy"  E: 


//_"  |x(T,j8)|2drd0  =  ff~  |u(t)|a|u(t +r)|adtdr  *  (2E)a  (4.64) 


Hence,  any  volume  moved  from  within  to  without  the  domain  of  physi¬ 
cally  interesting  values  of  r  and  (3  will  diminish  that  which  remains. 

It  turns  out  that  the  number  of  constants  C|j  (i,  j  =  1,  2)  can  be  re¬ 
duced,  so  as  to  simplify  the  algebra,  without  jeopardizing  the  result 
sought.  It  is  noted  first  that  the  cj0  enter  only  in  Eq.  4.58,  and  that 
these  constants  will  therefore  play  no  role  in  the  displacement  of  side 
lobes.  Stated  otherwise,  finite  discontinuities  in  the  phase  modulation 
4>(t)  are  irrelevant  to  the  displacement  of  side  lobes.  But  since  a  dis¬ 
continuity  in  4>(t)  is  nonphysical  anyway,  the  cj0  will  be  retained  for  the 
purpose  of  keeping  <}>(t)  formally  continuous.  The  Cjj  will  be  retained 
for  the  purpose  of  studying  the  effect  of  a  finite  discontinuity  in  4>'(t), 
the  instantaneous  frequency,  noting  that  such  discontinuities  seem  to 
be  in  keeping  with  fundamental  limitations  upon  the  signal.  Of  the  re¬ 
maining  constants,  either  the  ci2  or  the  ci3  can  be  put  equal  to  zero 
without  prejudicing  our  objective.  Under  the  second  of  these  alterna¬ 
tives,  the  simpler  theorems  of  Section  2  are  not  applicable  because 
924*n(t;  r,  P)/9t2  and  924>22(t;  t,  P)/9t2  vanish  identically,  and  the  side 
lobes  associated  with  4>j|(t;  t,  p)  and  4>22(t;  t,  P)  manifest  themselves 
only  as  straight  lines,  as  in  Eq.  4.7.  So,  for  the  present  we  put 

cia  =  csc  =  0  (4.65) 

The  alternative  is  considered  in  the  next  example. 

Equations  4.61  and  4.62  now  read 


94>1a(t;r ,  6)  _ 

at 


[(cll  ”  C21  “  C232’  ) 


2c23Tt  +  (c13  -  cS3)t?]  -  P  =  012(t ;t)  -  p  (4.66) 


23 


-  -2c23T  +  2(c13  -  c23)t 


(4.671 


d£«l>ip(t  ;t,/3) 

at2 


The  terminal  curves  associated  with  Eq.  4.56  are  given  by 

IT )  =  )]  t  =0  =  eu  -  csi  -  c23T‘  (4.68) 

j 

fit5 (T)  =  [®ia (t j ;t )]t.=T. t  =  c13(T  -  T)2  +  (ca  -  c2a  -  c23T2)  (4.69) 

with  T/2  <  t  <  T,  and  those  associated  with  Eq.  4.63  are  given  by 

gil)(T)=[0ia(tJ;T)]t.=T/2_T  =  cl3  (r  -  ~)  +  (cn  -  c21  -  cS3^)  (4.70) 

gJlV)  =  [®iS(tj;r )]t^mT/2  =  -c23(r  +  j)  +  (cu  -  c21  +  c13  —)  (4.71) 

with  0  <  t  <  T/2.  These  curves  are  sketched  in  Fig.  6  and  7  under  the 
additional  conditions 


Cfi3  <  C13  <  0 

(4.72) 

hi  +  (ci3  -  c33)  (— )  s  1/3  l» 

(4.73) 

where,  for  convenience,  we  have  put 

hi  =  cxi  -  cai 

(4.74) 

The  shaded  regions  (side  lobes)  are  determined  as  before  by  the  be¬ 
havior  of  [924>|2(t;  t,  P)/9t2]|=t.f  Eq.  4.67,  on  the  appropriate  t-  and 
T-intervals. 

The  inequality  (4.73)  requires  that  the  frequency  4>'(t)  be  discon¬ 
tinuous  at  time  t  =  T/2.  To  see  this,  suppose  to  the  contrary  that 
continuity  obtains,  ' ( T / 2 )  =  <j>2'(T/2).  Then 


.  whirh  case  Eq.  4.73  is  not  satisfied,  and  the  point  P  in  Fig.  6  coin¬ 
cides  with  the  origin.  If  the  required  jump  in  frequency  is  denote  y 

(A*)l2»  2 

nth,  =  *;(!)  =  -bl  ■ |C13'  C33’w) 

and,  upon  imposing  Eq.  4.73, 

(Sf)18  =  -|/3|.  <4'76) 

Having  seen  the  effect  of  continuity  in  frequency,  we  may  now  put 

cu  =  c21  =  0 

Upon  combining  previous  results,  we  have 

0x(t)  =  c10  +  ■“  Cist3.  0  s  t  *  ~l 


(4.77) 


0(t) 1 


with  Cjo  and  c20  chosen 
Cj3  and  C23  satisfying 


(4.78) 


-  fee-t3.  l<t«T 


0a(t)  =  cao  +  3  >  i 


to  obtain  continuity  in  phase  at  t  -  T/2,  and  with 


<  ci «  <  0 


C23  c13 


(4.79) 


(c13  -  c23)  (  2)  2 


(4.80) 


The  side  lobes  corresponding  to  the  first  and  third  integrals  in 
^  aLTo  to  «!>,  and  4>„,  are  also  of  interest.  By  previous  re- 

Eq.  4.55,  1.  e.  ,  11  ?  (4  791  they  appear  as  sketched  in 

Pigl  ^tundi^-rves  are  again  arcs  of  parabola,  The  fact 

thft  the  two  lobes  do  not  overlap  ^  the  region 

the  ambiguity  outside  the  central  lobe  of  IX(t,  PP  reduced, 

physically  interesting  values  of  x  and  p  the  central 

The  complete  ambiguity  Vagram  (fo  0,  now  equal 

lobe)  is  obtained  by  superposition  of  Fig.  b  tnrougn  o, 

to  zero. 
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Example  5. 


The  alternative  (and  possibly  more  practicable)  choice  of  the  con¬ 
stants  c2j»  c3j  discussed  just  prior  to  Eq.  4.65  is  now  discussed.  At 
the  same  time,  we  shall  subdivide  the  signal  duration  T  into  more 
parts  (four)  with  the  intention  of  displacing  more  of  the  ambiguity  vol¬ 
ume  beyond  ±|p|m,  and  of  distributing  in  azimuth  those  side  lobes  that 
are  connected  to  the  origin.  Thus,  we  have 


0(t)  = 


01  (t)  =  Cio  +  2  s  s  4 


02 (t)  -  Cso  +  2  C22ta>  4  S  i  2 


1  a  T 


3T 


03(t)  =  c30  +  2  c3at  >  2  S  ^  ^  4 


0*(t)  -  C4P  +  2  C4gt  ,  4  ^  t  s  T 


(4.81) 
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The  correctness  of  (4.87)  and  (4.88)  may  be  verified  by  means  of  a 
sketch  similar  to  Fig.  5. 

The  phase  functions  4>a  correspond  to  linear  FM  and  lead  to  the  de¬ 
generate  side  lobes. 


P  =  "ci2r 


(4.89) 


(see  Eq.  4.7  and  the  attendant  discussion.  )  It  turns  out  that  we  can 
choose  two  of  the  ca  to  be  negative  (and  unequal)  and  two  to  be  positive 
(and  unequal).  By  so  doing,  the  ambiguity  volume  near  the  origin  will 
tend  to  be  distributed  in  azimuth,  rather  than  concentrated  somewhere 
so  as  to  result  in  a  prominent  ridge  or  peak. 

The  side  lobes  corresponding  to  the  phase  functions  fy,  i  <  j,  can, 
at  the  same  time,  be  displaced  outward  beyond  ±|{3|m  .  The  terminal 
curves  are  given  by 

f| [V(T)  =  [«u(t;r)]tHr/4^  =  -c«T  +  (c12  -  Cj8) 


fif(T)  =  [®u(t;r)]tmlT/4  =  -C*r  +  (cia  -  Cj3)  ^ 


Under  the  conditions 


Cjl2  ^  0,  Coo  <  0,  C32  >  0,  C42  ^  0 


(cia  -  C32)  —  ^  |0|, 


(C48  "  c»)  ^  2  \P\w, 


(4.90) 


(4.91) 


(4.92) 


(4.93) 


Eq.  4.85,  4.89,  and  4.90  yield  the  partial  ambiguity  diagram  sketched 
in  Fig.  9-  The  dashed  lines  are  the  loci  of  Eq.  4.89;  the  dotted  lines 
represent  the  continuation  of  the  displaced  side  lobes  into  the  next 
interval,  T/4  <  t  <  T/2. 

The  remaining  T-intervals  involve  equations  similar  to  those  con¬ 
sidered  above.  Only  one  additional  condition. 


T  a 
C3B  ~2  S  I  P  I  * 


(4.94) 


it 


is  required  to  insure  that  the  remaining  side  lobes  are  outside  the 
range  of  physically  interesting  values  of  There  are  three,  they  all 
lie  below  —  |  p |m ,  and  each  has  the  character  of  the  lowermost  aid"  lobe 
shown  in  Fig.  9,  except  that  it  is  displaced  to  the  right  by  one  or  two 
multiples  of  T/4.  The  area  of  the  displaced  side  lobe  4>y,  i  <  j,  is 
I Cj2  -  cj2|(T/4)2. 

For  the  case  of  an  a(t)  which  is  nearly  rectangular  and  normalized 
to  unity,  A(tj;  t)  *  1/T,  and  the  displaced  ambiguity  volume  computed 
according  to  Eq.  3.17  under  the  assumption  that  the  side  lobes  are  dis¬ 
joint  is  easily  found  to  be  3/8  of  the  total  (unity).  If  the  number  N  of 
subdivisions  of  the  signal  duration  T  is  decreased  to  2,  the  displaced 
volume  is  1/4;  if  N  is  increased  to  8,  it  is  7/15;  and  if  N  is  increased 
indefinitely,  it  tends  toward  1/2.  It  turns  out  that  the  fraction  of  the 
displaced  side-lobe  area  which  lies  to  the  left  of  t  =  T/2  increases  as 
N  increases.  Hence,  an  A(tj;  t)  which  decreased  with  increasing  t 
would  result  in  a  greater  displaced  volume. 

The  question  of  whether  other  phase  modulations  will  yield  results 
which  are  more  favorable  in  the  sense  of  displaced  ambiguity  volume 
has  not  been  pursued. 

Example  6. 

Although  the  finite  jumps  in  the  instantaneous  frequency  <J>'(t)  which 
were  involved  in  the  previous  two  examples  do  not  seem  to  violate  any 
fundamental  restrictions  upon  the  signal,  they  are  likely  to  be  difficult 
to  realize  in  practice.  It  is  therefore  proper  to  ascertain  the  conse¬ 
quences  of  replacing  the  jumps  by  continuous,  rapid  transitions.  To 
this  end,  we  consider  the  modulations 


0.(t) 


0!  (t)  =  cie  +  |  c^t2,  Osts| 
03(t)  =  c*  +  |  c»t2,  |s  t  s  T 


(4.95) 


which  will  be  made  to  have  a  frequency  jump  at  t  =  T/2,  together  with 
the  modulation 
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0x(t)  =  c10  +  J  ciata, 

T 

0  s  t  -  c 

0b(t)  = 

^s(t)  =  C20  +  cait  +  J  C^t2, 

T  T 

7" €  st*I+c 

(4.96) 

tf^t)  =  c30  +  J  C$t2, 

for  which  the  instantaneous  frequency  will  be  made  continuous  by  vir¬ 
tue  of  a  linear  transition  during  the  interval,  T/2  -  t  <  t  <  T/2  +  c, 
with 


0<^«1  (4.97) 

The  constants  ci0  are,  again,  supposed  to  be  chosen  so  that  4>a(t)  and 
<j>b(t)  are  continuous.  In  the  case  of  4>b(t),  the  imposition  of  continuity 
in  frequency  yields 


ca  =  ■  h  (f 4 ')  (I ' c) ,Ca> ' Ca) 

=  -c»(f-c)]  <4-98) 

The  constants  c,2.  c32  are  chosen  so  that  the  degenerate  (linear  FM) 
side  lobes  corresponding  to  <J>u  and  $33  have  opposite  slopes,  and  the 
remaining  side  lobe  for  <J> a ( t)  is  displaced  beyond  +ipim: 

ci2  >  0,  C32  <  0 ,  -C33  ~  ^  1/3 1 a  (4.99) 

* 

For  definiteness,  we  also  take  |c?2|  <  |c32|.  The  resulting  ambiguity 
diagram  for  4>a(t)  is  sketched  in  Fig.  10;  as  before,  the  linear  FM  side 
lobes  are  indicated  by  the  dashed  lines. 

The  labeling,  fj^x)  and  g'.^x),  of  the  terminal  curves  employed  in 
this  example  means  that  the  associated  phase  function  is  4>jj,  with  i  <  j, 
and  that  this  is  the  k^  time  <|>. j  has  occurred  as  t  -*■  T  from  t  =  0  +  .  The 
functions  ^(t)  correspond  to  the  lower,  and  the  functions  g^(x)  to  the 
upper,  limits  of  the  integrals  with  which  they  are  connected. 
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FIG.  10.  Ambiguity  Diagram  for  the  Modulation  4t(t),  Eq.  4.95,  and  the  Conditions  (4.99). 

With  the  conditions  (4.99)  supplemented  byEq.4.98,  we  obtain  for  <|>b(t) 
the  side  lobes  sketched  in  Parts  a  and  b  of  Fig.  11.  Comparison  of  this 
figure  with  Fig.  10  shows  that  the  replacing  of  the  frequency  jump  at 
t  =  T/2  by  a  linear  transition  has  had  the  undesirable  result  of  generat¬ 
ing  additional  side  lobes  (the  linear  FM  lobe  4>22  in  Part  b  of  Fig.  11,  and 
both  of  the  side  lobes  in  the  two  parts  of  Fig.  11),  portions  of  which  lie 
within  the  domain  of  expected  values  of  t  and  (3.  However,  under 
Eq.  4.97,  and  to  within  the  SPA,  contributions  to  Ix(t,  P)|  of  these  ad¬ 
ditional  side  lobes  will  be  small  compared  to  lx(0,  0)  (  ,  as  will  now  be 
shown. 


For  the  side  lobes  <t>|2  and  4>23  we  have,  using  Eq.  4.98, 

83^T-,^=Cia-c.^(f4c)(c,a.c,) 

a  *"t=T'l£1  =  c,s  '  CsE  =  27  (I ' ') ,C!B " c,s) 


and  by  Eq.  3. 17, 


lxia(T,/3)|  *  — - JjP - —  A(t(3l2);r) 


,i/a 


/  T  \ 

\T +  c/  <Cls  ” c$)' 

(2f)1/s 

(f-f)  (ci2  -  ca)172 
where  the  stationary  points  are  given  by 


lxa,(T,/3)b 


A(t(23);r) 


f(l3)  - 


(i +  * )  (i " *  )(Cia  ~ Cse)  +  [Cae  ( I +  0  ~ Cia  (I  ~ c)] 

(j  +  «  )  (Cl2  "  C*) 


=  |  T  +°(f)'  °<TSf +£ 


and 


t(aa) 


(•J  +  e)  (y  -  c)  (c$  -  C1S)  +  (csT  +  /3)(2c) 

(~2  ~  *)  ^ C32  "  C12) 


=  T  +  0(t)’  0<TiT+c 


For  the  displaced  side  lobe  <J>13,  we  have 

d3<t>,  c(t ;T,  fl)  _ 


at£ 


-  cia  "  c33 


Ixiafr'.gll”  (;w  .!e^)175Alt)13>^') 


(4.100) 


(4.101) 


T  +  20C 


(4.102) 


(4.103) 


(4.104) 
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and 


t(i3)  =  csTL+P  2f  s/s  T 

cls  -  Cjb 


(4.105) 


From  Part  b  of  Fig.  11  we  find  (by  requiring  the  ordinate  p  to  lie  within 
the  shaded  region  4>13)  that 

\  -  t'  +  o(^)  s  tj13)  *  |  +  0  (^)  for  2f  sr'^+f  (4.106) 

and  that 

0  s  t(jl3)  5  T  -  r'  for  +  e  sT's  T  (4.107) 

(12) 

It  is  now  shown  that  if  t  and  tj  are  prescribed  according  to  Eq.  4.102, 
then  a  c-13)  and  a  t'  satisfying  Eq.  4.106  can  be  found  such  that 

A(t(4ls);r)  =A(t(J13);r/)  (4.108) 

We  use  the  continuity  of  a(t)  and  work  only  to  within  terms  which  are 
0(«/T).  By  doing  the  latter,  small  areas  of  the  side  lobes  are  excluded; 
these  can  be  taken  into  account  by  a  continuity  argument.  We  have 

A(t$ls)  ;t  )  =  a(tSl8,)a(t(Jw)  +  T)  «  a  -  r)  a  f|)  ,  0  <  r  <  \ 

Since  a(t)  is  continuous  and  vanishes  for  t  >  T,  there  exists  a  t, 

0  <  t  <  T/2,  such  that 

=  a(^  +  f)  a  f|)  ,  0  <  f  <  j  (4.109) 

Now,  Eq.  4.106  is  satisfied  by  tj13)  =  T/ 2  and  t'  =  t,  and  for  these 
values  we  obtain  Eq.  4.108.  The  function  A(tj23)  ;  t),  where  tj23)  and  t 
satisfy  Eq.  4.103,  is  given  by 

A(tf3);T)  =  a(tjS3) )a(tj33)  +  T )  »  a  (j)  a  (j  +  r)  ,  0  <  r  <  \ 
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Upon  comparing  this  with  Eq.  4.109,  it  is  seen  that 


A{tf3);T)  =A(t 


(13)  .T'v 
)  >T  ) 


(4.110) 


for  tj13)  -  T/2  and  t*  =  t.  Next,  the  approximations  in  Eq.  4.108  and 
4.110  are  used  in  Eq.  4.100  and  4.101,  thus  arriving  at 


lxjr.fi)!  *  2(|)‘'2|x,3(r'  fi')\ 


by  way  of  Eq.  4.104.  Finally,  we  employ  the  general  and  exact  rela¬ 
tionship, 

lx(o,  0)| 

(4.111) 

to  obtain 

|XIS(T,/9)|  <  2  (J 

-V  1/2 

lx(0, 0)1 

K»(T’^)i  s  2  (t 

\l/2 

-)  lx(0,0)| 

(4.112) 

as  claimed. 

The  demonstration  that  Ix22(t,  P)|«Ix(0,  0)|  follows  along  similar 
lines.  Application  of  Erdelyi's  theorem  to 


li 


(T',  P')  =  f  /2tT  A(t;T/)e2’,i*1,(t:T',',0dt 
Jo 


and 


gives 


Jr  12-f 


'll 


2jr|c13r'  +  p'\ 


h,y,m\  - ♦  o[^] 
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* 


and 


a2(I+ f  • TiT)  +A’( * ' f :T) 


2a(^  +  €  -  t;t)a(:|  -  €;t)  cos  2ff(caaT  +  0)(2c  -  T) 


where 


and 


i3'  ^  ci3t/ 


0  <  r  <  2e,  /3  i  caaT 


°[{Ewcaa)2] 

(4.113) 


(4.114) 


Upon  using  Eq.  4,97.  4.114,  and  the  second  of  Eq.  4.98.  we  obtain 

4,’(|) 

lx.(T.|J)|-f- - — 5— |.in»(c«r+ff)P«  -T)l 

ff  |  (c33  "  Cja)T  +  rr,  P  I 


IxjT.fl)!  5  T 


4a‘ 


(1) 


and 


jXaa 


Xa(T,glLs  ± 

X^.fjT  T 


ff  1  (C3J  ~  Cia)T  +  rj.  $1 

8a  (-|)  Icxal'  +  0'l 
l  a  -  t')  |  (C33  “  cia)T  +  ^0l 


Given  t  and  (3,  it  is  always  possible  to  find  a  r'  and  a  p'  ******** 
(4, 113)  such  that  the  quantity  within  the  braces  is  unity.  Hence,  by 


means  of  Eq.  4.111,  we  arrive  at 

.  flv|  £  Jxaailill 
'x2a(T’/3)1  iv  (T'./n 


j\(0,  0)|  ^  3T  |x(°*0)i  (4.115) 


as  was  to  be  shown. 

Phase  modulations  of  the  sort  considered  in  this  anc I  the 
example  have  also  been  discussed,  in  a  different  way,  by  Rihacaek  an 

Mitchell  (Ref.  18). 


39 


5.  SUMMARY  AND  DISCUSSION  OF  RESULTS 


It  has  been  shown  that  the  stationary  phase  approximation  (SPA)  of 
Woodward's  time-frequency  correlation  function  x(t,  3),  Eq.1.1,  can 
be  elaborated  upon  so  as  to  yield  inequalities  (or,  alternatively,  fami¬ 
lies  of  curves)  specifying  those  domains  of  the  time-frequency  (t  -  (3) 
plane  on  which  the  ambiguity  function  Ix(t,  (3)|2  is  relatively  large 
(side  lobes).  The  elaboration  consisted  of  working  a  condition  satis¬ 
fied  by  the  stationary  points  for  large  Ix(t,  (3) I2  into  one  satisfied  by 
(3  and  t.  A  brief,  general  discussion  was  given  in  Section  3,  and  some 
examples  were  worked  in  Section  4.  The  method  fails  in  some  small 
neighborhood  of  the  origin  (r  =  0  =  (3),  but  the  behavior  of  Ix(t,  p)|2 
there  can,  in  principle,  be  determined  by  means  of  a  Taylor  series 
expansion.  Although  the  SPA  theorems  which  are  usually  used  do  not 
apply  along  tne  p-axis,  a  more  general  theorem  given  by  Erdelyi 
(Ref.  16  and  Section  2)  can  be  used  to  study  the  ambiguity  function 
there;  alternatively,  Eq.  3.24  can  be  used.  For  modulations  such  that 
the  amplitude  modulation  a(t)  varies  slowly  compared  to  the  phase 
modulation  <j>(t)— except,  possibly,  at  the  beginning  and  end  of  the 
pulse  — this  technique  should  be  a  useful  tool  for  signal  design. 

In  the  examples,  we  considered  first  those  phase  modulations  which 
satisfy  the  hypothesis  of  the  SPA  theorems  usually  employed  (Ref.  15 
and  Section  2):  i.  e.  ,  functions  <j>(t)  which  are  analytic  over  the  duration 
T  of  the  signal.  The  same  was  assumed  concerning  the  amplitude 
modulation  a(t).  The  weaker  hypothesis  of  Erdelyi's  theorem  was  thus 
satisfied  at  the  same  time.  It  was  found  that  in  these  cases  there  were 
two  side  lobes  (one  obtainable  from  the  other  by  inversion  in  the  origin) 
connected  to  the  origin  and  extending  outward  to  t  =  ±T. 

We  then  considered  phase  modulations  <p(t)  in  the  form  of  poly¬ 
nomials  which  satisfied  the  theorems  only  piecewise,  with  a(t)  the 
same  as  before,  and  it  was  found  that  (1)  some  of  the  side  lobes  could 
be  displaced  outward  along  the  p-axis  beyond  the  largest  expected  val¬ 
ues  of  p,  ±|ptm,  and  that  (2)  the  remaining  side  lobes  could  be  distrib¬ 
uted  in  azimuth  about  the  origin.  In  connection  with  Item  1,  finite 
discontinuities  in  <p(t)  were  found  to  be  irrelevant,  and  finite  discon¬ 
tinuities  in  4>' (t),  the  instantaneous  frequency,  were  found  to  be  essen¬ 
tial.  Because  of  the  volume  invariance  property  of  Ix(t,  P)|2,  Eq.  4.64, 
signals  of  this  sort  raise  the  hope  of  approximating  a  qualitatively  ideal 
ambiguity  function,  i.  e.  ,  one  having  a  sharply  peaked  central  lobe  and 
very  low  values  over  the  remainder  of  the  rectangle  bounding  the 


expected  values  of  t  and  p.  Although  the  existence  of  finite  jumps  in 
the  instantaneous  frequency  <j>'(t)  does  not  seem  to  stand  in  contradic¬ 
tion  to  fundamental  reJtrictions  upon  the  signal,  it  is  possible  that 
they  are  difficult  to  realize  in  practice.  For  this  reason,  we  con¬ 
sidered  in  the  last  example  a  phase  modulation  such  that  the  jump  was 
replaced  by  a.  linear  transition  occurring  in  an  interval  c.  We  found 
that  the  contributions  to  Ix(t,  (3) I2  of  the  additional  side  lobes  generated 
by  the  transition  were  small  compared  to  lx(0,  0) I2,  provided  that 
«  «  T. 


In  cases  in  which  a  large  portion  of  the  ambiguity  volume  corre¬ 
sponding  to  a  given  signal  already  lies  beyond  ±lp|m,  not  much  will  be 
gained  by  the  techniques  of  Examples  4,  5,  and  6  unless  the  signal 
duration  T  is  divided  into  many  (equal)  parts.  The  result  of  division 
will  be  a  signal  which  is,  perhaps,  undesirably  complicated,  or  which 
has  too  large  a  bandwidth.  Let  us  recall  that  Eq.  1.1  rests  upon  the 
narrow- band  approximation,  and  that  this  approximation  may  be  a 
poor  one  for  modulations  such  that  the  SPA  is  valid. 

Finally,  it  is  pointed  out  that  specific  modulations  that  are  arrived 
at  by  means  of  the  methods  described  here  are  best  considered  as 
candidate  signals,  because  we  have  made  use,  both  explicitly  and  im¬ 
plicitly,  of  just  the  leading  term  of  an  asymptotic  expansion.  While 
asymptotic  expansions  can  yield  very  accurate  results,  it  is  necessary 
to  verify  their  accuracy  by  looking  at  more  terms  in  Eq.  2.7  and  2.8. 
Alternatively,  one  can  resort  to  numerical  integration  of  Eq.  1.1,  or, 
possibly,  to  analog  simulation  of  the  matched  filter  which  is  charac¬ 
terized  by  x(t,  p). 
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